Temperature dependence of uniform static magnetic susceptibility in a 
two-dimensional quantum Heisenberg antiferromagnetic model 
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A perturbation spin-wave theory for the quantum Heisenberg antiferromagnets on a square lattice 
is proposed to calculate the uniform static magnetic susceptibility at finite temperatures, where a 
divergence in the previous theories due to an artificial phase transition has been removed. To 
the zeroth order, the main features of the uniform static susceptibility are produced: a linear 
temperature dependence at low temperatures and a smooth crossover in the intermediate range and 
the Curie law at high temperatures. When the leading corrections from the spin-wave interactions 
are included, the resulting spin susceptibility in the full temperature range is in agreement with the 
numerical quantum Monte Carlo simulations and high-temperature series expansions. 
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I. INTRODUCTION 

The quantum Heisenberg antiferromagnet (QHAFM) 
is a prototype model to describe the magnetic proper- 
ties in the parent compounds of high temperature super- 
conducting cupratesi Non-linear sigma model is shown 
that its ground state has long-range order and its low en- 
ergy excitations are in a renormalized classical state^ 3 -^ 
These results are confirmed by the quantum Monte Carlo 
methods^ and are consistent with experimental data. 2 
Moreover, the Schwinger boson mean field theory^ and 
the variational spin-wave theory* have also been devel- 
oped for the two-dimensional QHAFM, leading to sim- 
ilar conclusions. Previous theoretical studies have con- 
centrated in the low temperature region. However, the 
low-energy interesting physics of the two-dimensional 
QHAFM can persist up to the temperature region T ~ J. 
Moreover, the anomalous normal state in the charged 
doped QHAFM in the cuprates^ relies on a theory for 
the magnetic fluctuations in the intermediate tempera- 
ture range. Therefore, a theory is needed to be well de- 
fined in the full temperature range. 

Although the Schwinger boson mean field theory^ 
and Takahashi's variational spin-wave theory^ for a 
QHAFM in square lattices can be extended to the fi- 
nite temperature range, the mean field ansatz results 
in a finite temperature phase transition at T c — 0.91J, 
which explicitly violates the Mermin- Wagner theorem 11 . 
To overcome this artificial mean filed divergence of the 
magnetic susceptibility, we will propose a perturbation 
spin-wave theory. 

In the Takahashi's variational spin-wave theory, the 
Dyson-Maleev spin-boson representation is used for the 
QHAFM on a square lattice with a Neel ordered ground 
state: 



S t = a], = (25 - a\a,i)ai, S? = S - 
for the spins in the spin-up sublattice, and 



a\ai, 



Sr = b h S+ = b](2S - b]bj),S? = -S + b]b 3 , 



(1) 
(2) 



for the spin-down sublattice. The Dyson-Maleev or 
the Holstein-Primakoff spin-wave theory works well for 
the ground state. The corresponding results, such 
as the spin-wave velocity, the spin stiffness constant, 
the sublattice magnetization, and the perpendicular 
susceptibilit y 12 ' 13 ! 14 agree well to the quantum Monte 
Carlo simulations.— 

To develop a finite temperature spin-wave theory 
Takahashi introduced a constraint for the vanishing local 
magnetization, i.e., 



(3) 



at the i'th lattice site, which automatically fulfills the 
Mermin- Wagner theorem. Takahashi also introduced a 
mean field order parameter {a%bj) 



b\) to account for 



* 3 

the spin-wave interaction in mean field level. This mean 
field ansatz can describe reliably the low temperature 
physics, but fails at high temperatures, where it produces 
an artificial phase transition. 

In this paper, we present a perturbation spin-wave 
theory with the constraint Eq. ([3]). The spin- wave in- 
teractions are considered within a many-body perturba- 
tion method. The mean field divergence in the previous 
variational spin-wave theory is removed. To verify our 
perturbation spin-wave theory, we calculate the uniform 
static susceptibility. Its temperature dependence in the 
zeroth order can catch most features given by the quan- 
tum Monte Carlo simulations and high temperature se- 
ries expansion. It shows a linear temperature dependence 
at low temperatures, a smooth crossover in the interme- 
diate range and the Curie law at high temperatures. By 
including the first order corrections from the spin wave 
interactions, our perturbation spin-wave theory agrees 
quantitatively well with the numerical simulations. 

This paper is organized as follows. In Sec. [Til we 
present the linearized spin wave results up to the zeroth 
order. In Sec. IIIH a perturbation theory with the first or- 
der corrections from the spin-wave interactions are devel- 
oped. We also give the comparison of the static magnetic 
susceptibility with the results of quantum Monte Carlo 
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simulations and high-temperature series expansions. In 
Sec. IIV[ we give the discussion and summary. 



II. LINEARIZED SPIN- WAVE THEORY 

In this paper we mainly discuss a QHAFM on a quare 
lattice with Hamiltonian defined by 



H = J^S 4 



(4) 



where J > and (ij) denotes nearest-neighbor sites. 
Since the ground state of this QHAFM Hamiltonian is 
an Neel ordered state, we separate the square lattice with 
2N sites into sublattice A for the up spins and sublattice 
B for the down spins. 

The Dyson-Maleev spin representation given by Eq. |T]) 
and (J2J) are defined upon the Neel ordered state. In terms 
of boson operators, the model can be divided into 



H = H l +H I , 



(5) 



where Hi contains the classical energy and the quadratic 
terms, 



Hi = 2JSY^ (<4 a i + b ] h 3 + a i b j + 4 b j) 

(ij) 

-NJzS 2 , 



(6) 



and Hi contains the quartic terms, describing the spin- 
wave interactions, 



H I = -jJ24(a l + b]) b, 

(ij) 



(7) 



In this section, we neglect the spin-wave interactions 
and focus on the Hamiltonian Hi . We refer this approxi- 
mate to a linearized spin- wave theory (LSW). The spin- 
wave interaction in the Hamiltonian Hi will be studied 
in Sec. ED 

The constraint Eq. ([3]) , from the Dyson-Maleev repre- 
sentation, can be described by introducing a Lagrangian 
Hamiltonian as^ 



H> 



E 



Ui(S - a\ai) 



3 eB 



where the Lagrangian multipliers fii and fij are assumed 
to be site independent, i.e. /i^ = /i 3 = /i. Physically 
the Lagrangian multiplier /i plays a role of an effective 
chemical potential for the spin waves and leads to a finite 
gap, thus an effective mass, for these bosonic excitations 
at finite temperatures. For convenience in our following 
discussion, we define /i = JzS(\ — 1) with the coordinate 
number z = 4. 

By using the Fourier transformations and performing 
the Bogoliubov transformations, the Hamiltonian Hq = 



Hi + H\ can be expressed as 



H = £ k (4"k + 4Ac) + E £k 

k k 

-NJz\S{l + 25) + NJzS 2 , 



(9) 



where £k = JzS^J A 2 — 7^ and 7k = cos \k x cos \k y . 
Here ^aj^ak) and (^P^Pk^j are operators describing two 
branches of Bogoliubov excitations. The free energy per 
site is given by 



k 

JzS 2 



(10) 



Minimizing the free energy with respect to A leads to a 
self-consistent equation 



s+ 9 - ?^fE 



A 



coth 



2 2N ^ y/\z _ j£ ' \2TJ 



(11) 



as a consequence of the constraint Eq. (|3|). 

The uniform static susceptibility per site can be calcu- 
lated from the static spin-spin correlation function, 



+r)l- 



(12) 



Here the thermal average is defined in the zeroth 
order linearized spin- wave approximation as (O)i = 



It can be easily shown that the transverse 



Tr(e- H o/ T 6) 
Tr(e- H o/ T ) 

spin-spin correlation function x_l = J2WT Y^iA^l ^l+r + 
SfS^ r ) = 0. Therefore, only the longitudinal spin- 
spin correlation function Xz = o^T^2i r (^i^i+r) nas a 
finite value. The vanishing transverse spin-spin correla- 
tion function can also be shown in the Takahashi's vari- 
ational spin-wave theory^ Physically it stems from the 
local constraint given in Eq. (J3j> . 

The longitudinal spin-spin correlation function can be 
calculated by introducing a weak external magnetic field 
Be z along z-axis which introduces an additional term 

h b = -B (j2 ieA Sf + J2 j€ b Sj) t0 the model Hamil " 
tonian. The longitudinal static susceptibility thus reads 

Xz = ~ d Jb? ^ ; where the free energy f(B) includes the 
contribution from the external field Hb- Then the uni- 
form static susceptibility within the linearized spin- wave 
theory has the form: 



Xi(T) 



1 



12NT 



(13) 



A similar expression has also been obtained by Takahashi 
in his variational spin- wave theory^ 
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1. Antiferromagnetic ordered state at zero temperature 

The ground state of the QHAFM on a square lat- 
tice has a long-range order— ^ due to the presence of 
Bose-Einstein condensation (BEC) at the momentum 
k = (0, 0). Since A = 1 + O (N' 1 ) at zero temperature, 
the spin-wave spectrum is given by 

e k = JzSyJl-il. (14) 

Comparing to the Takahashi's spin wave spectrum- 
EifT,k = Jzm\y/\ — 7^ with mi — S + 0.078974, we 
find that it can be obtained by rescaling our spin-wave 
result as 

£MT,k = a£k, a = —. (15) 

The additional factor a is the main difference at low tem- 
peratures between our linearized spin-wave theory and 
Takahashi's variational theory . 

The spin-spin correlation function in the long distant 
limit can be derived as 

(S r S l+r ) = c(r) (m + , (16) 

where c(r) = 1 if I, l+r are in the same sublattice, and —1 
if I, l + r belong to different sublattices. mo = S — 0.19660 
is the magnetization at zero temperature, which is the 
same as that in the Schwinger boson mean field theory^. 
Both the spin-spin correlation function and the local 
spontaneous magnetization mo in our linearized spin- 
wave theory agree with those of the variational spin- wave 
theory^ These results imply that, to the order of 0(S^ 1 ), 
the spin-wave interactions have no contributions to the 
spontaneous magnetization and static spin-spin correla- 
tion at T — OK, consistent to the previous resultsj 12 ' 14 ! 17 
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FIG. 1: The Lagrangian multiplier A(T) in the linearized spin- 
wave theory for the QHAFM on a square lattice with different 
spin magnitude S = 1/2, 1 and 3/2. 
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T/JS(S+1) 

FIG. 2: The uniform static magnetic susceptibility \i{T) for 
S — 1/2, 1 and 3/2 in our linearized spin wave theory of the 
QHAFM on a square lattice. 



2. Paramagnetic state at finite temperatures 

At finite temperatures, according to the Mermin- 
Wagner theorem, the long-range order is destroyed in the 
QHAFM model by quantum fluctuations*^ Fig. Q] shows 
the temperature dependence of the Lagrangian multiplier 
A(T) with different quantum spins S = 1/2,1,3/2. At 
low temperatures T <C JzS, we follow the Takahashi's 
method and show analytically that A approaches expo- 
nentially to its BEC value 1 as 



exp 



irJzSmQ 
f 



(17) 



The uniform static spin susceptibility in this low temper- 
ature region can also be derived analytically, leading to 
a linear-temperature dependence, 



m 
3JzS 



2T 



3ir (JzS) 



(18) 



The first term gives a finite spin susceptibility at zero 
temperature. The linear term results from the linear den- 
sity of states of the spin-wave excitations in the low en- 
ergy limit. The low-energy spin- wave excitations at low 
temperatures have an approximate dispersion 



e k = JzS\ A 2 



4 k ' 



(19) 



where a dimensionless bosonic gap is defined by A = 
At low temperatures, A decays exponen- 
therefore the energy spectrum 
The density 



Va^T 

tially to 1 as Eq.(fT7] 
Eq. (fT9")) is approximately as £k = ^-k 
of states p(oj) of the spin-wave excitations is p(to) 



J d 2 kS (u 



It is this linear density 



Tt(JzSY 

of states that leads to the linear-temperature dependence 
of the static magnetic susceptibility at low temperatures. 
Note that by scaling the linear temperature susceptibility 
Xi(T) by the factor a, we can obtain the low-temperature 
susceptibility given by Takahashi^ xmt(T) = \xi )• 
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At low temperatures, the spin-spin correlation function 
follows 



(Si • S i+r ) = c(r) 



T 



7r \ JzS J r 



/€ 



(20) 



where c(r) is defined in Eq. (fT6]) and £ is the correlation 
length defined by £ = 4v / A \_ 1 ■ From Eq. (fTT|) , it can be 
easily shown that 



, f rr\ JzS ( ^JzSniQ 
£( T ) = exp 



4T 



2T 



(21) 



Takahashi's result for the spin-spin correlation length 
£mt (T) can be obtained by rescaling with the factor a de- 
fined in Eq. (fl5|) . !;mt{T) = £ (^). This exponential tem- 
perature dependence of the spin correlation length agrees 
with that obtained from one-loop renormalization group 
approach of the non-linear sigma modeli^S Compared 
to the experimental data and the numerical simulations, 
the two-loop approximation is necessary to convert the 
1/T prefactor of the exponential into a constant!^ 

When increasing the temperature, the Lagrangian mul- 
tiplier A(T) in Fig. [Devolves from the low-temperature 
exponential behavior to the linear temperature depen- 
dence at high temperatures. At the same time, the uni- 
form susceptibility as shown in Fig. [2] increases firstly 
to a broad peak at the intermediate temperature T c oc 
JS(S + 1), and then decreases at high temperatures. 

At high temperatures T 3> JzS, the dispersion of the 
energy spectrum becomes relatively weak and A > 1. 
The self-consistent equation Eq. fTTj) gives rise to 



A(T) 



T 
7z~S 



In 1 



(22) 



In the high-temperature region, the uniform static mag- 
netic susceptibility obeys the Curie law, 



Xi(T) = 



(S + 1)S 
3T ' 



(23) 



Such a high-temperature Curie-law behavior in the lin- 
earized spin-wave theory is exactly the same as that in 
the Takahashi's theory^ Why our linearized spin-wave 
theory can give rise to a Curie-law behavior, in contrast 
to the general spin- wave theory is a very interesting ques- 
tion. The basic reason is that there is an effective chem- 
ical potential introduced for our spin waves by the local 
constraint Eq.f3]). This gives rise to a temperature de- 
pendent finite gap for the spin waves which behave as 
bosons with effective mass at finite temperatures. 

It should be noted that in our numerical solution to 
the self-consistent equation Eq. (|lip and in calculating 
the uniform magnetic susceptibility Eq. (|13[) . there is a 
trick to include the contribution of the low-energy spin- 
wave excitations at low temperatures. Although no exact 
BEC occurs at finite temperatures, in the actual numer- 
ical calculation Eq.fTTj) shows us that the exponential 
decay of the bosonic gap JzS A at k = (0, 0) will lead 
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FIG. 3: Comparison of the uniform static magnetic suscep- 
tibility with the 3 = 1/2 linearized spin- wave theory (LSW) 
and the Takahashi's variational theory 9 for the QHAFM on 
a square lattice. 



to a difficulty in numerically calculating the momentum 
integral from the low-energy region <C T. We call this 
difficulty as a BEC-like pseudo-singularity. 

A trick to deal with the BEC-like pseudo-singularity is 
to separate the momentum space into two regions: Re- 
gion (i) k £ k-BEC + <>k and Region (ii) is the rest part 
in the first Brillouin Zone. Here the momentum k € 5^ 
satisfies |k| < k m . k m is an irrelevant parameter and we 
choose k m = 0.01 in our calculations. The momentum 
integral in Region (i) can be firstly calculated with an 
approximate spin- wave energy dispersion Eq. (|19p . 

For comparison, the uniform static magnetic suscepti- 
bility from our linearized spin- wave theory and the Taka- 
hashi's variational results are shown in Fig. [3J The sharp 
peak value at T c = 0.91J in the Takahashi's theory comes 
from an artificial phase transition due to vanishing of the 
mean-field ansatz. At low temperatures, the static mag- 
netic susceptibility from our linearized spin-wave theory 
is slightly larger than the Takahashi's result. At high 
temperatures, the Curie law occurs exactly in both the- 
ories. Clearly, the artificial phase transition has been 
removed in the linearized spin-wave theory. 

One special feature in the uniform magnetic suscepti- 
bility shown in Fig. [2] is a kink structure around the tem- 
perature T w 0.5 JS, where the slop of the susceptibility 
changes from a smaller value into a larger one. This spe- 
cial feature is also observed in the quantum Monte Carlo 
dataJ^ It was argued to originate from a crossover from 
a renormalizcd classical regime into a quantum critical 
regime. However, the Monte Carlo simulation shows that 
a similar abrupt change also appears in a classical Heisen- 
berg model Therefore, whether there is a crossover is 
still controversial. However, this kink structure implies 
obviously different scenarios occurring in these two tem- 
perature regions. 
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III. FIRST ORDER CORRECTIONS FROM THE 
SPIN- WAVE INTERACTIONS 

In this section, we will go beyond the linearized ap- 
proximation by including the first order corrections from 
the spin-wave interactions of the QHAFM model. The 
method is based on the many-body perturbation theory. 



A. Formulation for first order corrections 

After the Fourier and Bogoliubov transformations, the 
total Hamiltonian H = Hi + H\ + Hj can be approxi- 
mated by2£ 

TL = Ti-o + Hi + E c , 
U a = ^2 A * ( nk + n k) ' 

k 

ki,k 2 

NJz 



E c = J2 £ k + NJzS [S - A (1 + 2S)] 



4, (24) 



k 

t. 



where n k = a k a k , n k = /3 k /3 k , ci(A) = 1 - £] k 



JzSCu ( 



l + Ci 



A — 



~k 



rep- 



Ck = \/A 2 - 7 k , and A k = 

resents the spin-wave energy spectrum with corrections 
from the normal ordering of Hi. At zero temperature, 
A = 1 + O (N^ 1 ) and the spin- wave energy dispersion is 

given by A k = JzSCk (l + ^) with c x (1) = 0.15794 

for S = 1/2. The parameters and B^> are defined 
by 



S (2) 




Cki Ck 2 

(A -7^) (A 



7k 2 



CkXi 



k 2 



In the above approximation, we have ignored the off- 
diagonal parts in the Fock space and picked up all the 
terms of the two spin- wave interactions. This is enough 
when we only consider the first order corrections from 
the spin-wave interactions for the free energy, since the 
linked cluster theorem tells that f = fo + (Hi)o, where 

Tr(e- H o' T d) 
Tr(e-"o/T) 

and fo is defined for the free bosons with energy spectrum 
A k . 

The free energy per lattice site reads 



the thermal average is defined by (O)o 



/ 



T 
N 



i 



-A k /T 



Jz 



1 



-.E r 



(25) 



where c 2 (A) = £ £ k ^T™k and ,~i k = (e A *' T - l) \ 
Minimizing this free energy, we can obtain the following 



self-consistent equation for the Lagrangian multiplier A, 



S + l = JLVoothf^ fA^(A-l) 7 2 
) o.N \ 2T J VCk 

(A - it) 



2T J V Ck ' 2SC£ 



1 x 



A z ^4Tsinh 2 (#)^Ck 9A ' 



(26) 



k " V22 

where c(A) = ci(A) — 2c2(A) and 



<9A 



^(A + C1 (A-1)7: 



-J; 



Ck ' 25^ 
A-7g^ (A-l)7g 2 
27VCk 



k 2 



It should be noted that the self-consistent equation 
Eq. ([26|) is physically equivalent to the first order per- 
turbation expansion for the constraint Eq.Q, 

= S-(al ai ) + [ 1 dr(T T W J (r)at(0)a i (0)) . (27) 
Jo 

The uniform static susceptibility can also be calculated 
from the spin-spin correlation function given in Eq. (|12j) . 
When the first order corrections from the spin- wave inter- 
actions are included, the transverse static susceptibility 
has a small but finite contribution, 



X±(T) = fX %(T), 



(28) 



XF (T) 



k 



c 3 



Ecoth 



2A 



A*\ fX , Cl (A-l) 7 2 



277 VCk 



25C, 



In the derivation for we have used the perturbation 
constraint equation Eq. (|27p and ignored the higher order 
terms of the spin- wave interactions. 

The longitudinal static susceptibility can also be calcu- 
lated from the free energy by introducing a weak external 
field Hamiltonian Hb, yielding 



X, (T) 



WfT,^ (#)+M^-c 2 ], (29) 



where the second term comes from the two spin-wave 
interaction compared to the linearized spin-wave expres- 
sion Eq. (fT3"|) , and 



c 3 (A) 



N 4-^ 



(A -ll) 



4TCkSinh 2 (^) 



1 V (A- 7 g)coth(#) 

- N 2^ 4T2CkSinh2( ^) ■ 

The uniform static susceptibility is thus given by 
X{T) = \[XX (T) + Xz (T)}. 



(30) 
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FIG. 4: The Lagrangian multiplier in our perturbation spin- 
wave theory with the first order corrections from the spin- 
wave interactions. 



B. Numerical results 

The self-consistent equation Eg. ((26)) and the static 
magnetic susceptibility Eq. (|30p have been numerically 
calculated. When the first order corrections of the spin- 
wave interactions are included, the temperature depen- 
dence of both the Lagrangian multiplier A and the uni- 
form static magnetic susceptibility have similar behaviors 
to the linearized spin-wave theory. 

The Lagrangian multiplier A is shown in Fig. [4] for 
different spin S = 1/2,1 and 3/2. At low temper- 
atures, it decays exponentially to its zero-temperature 
value 1, indicating the gapless spin wave excitations. At 
high temperatures, A has a linear-temperature variation 
A(T) = Ko+KiT, where kq and K\ have weak temperature 
dependence. For S = 1/2, kq = 0.334 and k\ = 0.583 in 
the temperature region 3 J < T < 10 J. 

The static magnetic susceptibility is displayed in Fig. 
[SJ It has three main features: (i) the linear-temperature 
behavior at low temperatures, (ii) a broad peak at an in- 
termediate characteristic temperature T c cx JS(S + 1), 
(iii) the Curie- Weiss law at high temperatures. For 
S = 1/2, a low-temperature linear behavior can be fitted 
by x(T) = Xc + bT, with X c = 0.0418 and b = 0.0471. 
Compared to the Takahashi's theory where a similar lin- 
ear behavior with Xc = 0.0437 and b = 0.0396,— our zero 
temperature susceptibility Xc is slightly smaller than that 
of Takahashi's. Also the slope b is also slightly larger than 
that of the Takahashi's. 

At high temperatures, the spin-wave interactions mod- 
ify the exact Curie law of the static magnetic susceptibil- 
ity Ea.(|23p into a Curie- Weiss behavior. As shown in Fig. 
El it evolves into the behavior obtained from the numer- 
ical high-temperature series expansions. Physical origin 
of the high-temperature Curie- Weiss law also comes from 
the linear-temperature dependence of Lagrangian multi- 
plier A in the energy spectrum of the spin waves. 

It should be noted that the kink structure in the static 




T/JS(S+1) 



FIG. 5: The static magnetic susceptibility with different spin 
S — 1/2, 1 and 3/2 in our perturbation spin- wave theory with 
the first order corrections from the spin-wave interactions. 




T/J 



FIG. 6: Comparison of the uniform static susceptibilities from 
the linearized spin-wave theory (LSW) and the perturbation 
spin-wave theory (PSW) with first order correction for the 
QHAFM on a square lattice. We also show the result from the 
Takahashi's modified spin-wave theory and numerical data 
from Quantum Monte Carlo and high-temperature series ex- 
pansions (HTSE). The inset shows the comparison of PSW 
with HTSE in high temperature region. QMC a comes from 
Ref. H and QMC 6 comes from Ref. 0. HTSE comes from 
Eq.(34) of Ref. & For comparison, QMC b have been multi- 
plied by a factor ~. 



magnetic susceptibility in the linearized spin-wave the- 
ory shown in Fig. [2] also occurs in the perturbation the- 
ory. This implies that the kink feature is a characteristic 
property in the QHAFM model. 

In Fig. [6j we compare the static magnetic susceptibil- 
ity of our spin- wave theory to the quantum Monte Carlo 
data£*^ and the high-temperature series expansions. It 
shows that our linearized spin- wave theory can reproduce 
all the main features in the numerical simulations: the 
low temperature linear behavior, the smooth crossover in 
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the intermediate temperature region and the Curie- Weiss 
law behavior at high temperature. This good agreement 
encourages us to study the frustrated J\ — J2 model pro- 
posed to describe the magnetic fluctuation in the newly 
discovered iron-pnictide superconductors i^i The first or- 
der correction from the spin-wave interactions makes the 
uniform static susceptibility both at low and high tem- 
peratures agrees quantitatively with the numerical quan- 
tum Monte Carlo and high temperature series expan- 
sions. This good fitting in Fig. [6] indicates that the 
physics of QHAFM model, at least to the static mag- 
netic response, has been captured by our perturbation 
spin-wave theory in the full temperature region. 

IV. DISCUSSIONS AND CONCLUSIONS 

The good agreement of our perturbation spin-wave 
theory with the numerical simulations shows the validity 
of our theory at finite temperatures. It seems surprising 
because the spin- wave theory developed in the pioneering 
paper o 22 ' 23 ^ 24 ' 25 is generally thought to work well only at 
low temperatures in the ordered phase. The main dif- 
ference of our finite temperature spin- wave theory is the 
requirement of the vanishing local magnetization. It is 
this local constraint that restricts the number of the spin 
wave excitations which thus behaves as bosons with an 
effective mass at finite temperatures. 

The extension of the spin-wave theory to the low- 
dimensional spin systems have been shown to work well. 
In one-dimensional quantum magnets, no magnetic order 
can survive at zero temperature due to the strong quan- 
tum fluctuations. Although there is no long range order 
in the ground state, the free energy and the static mag- 
netic susceptibility of the one-dimensional fcrromagnet 
from a modified spin-wave theory show same behaviors to 
Bethe-Ansatz solutions. 26 The Haldane gap for the one- 



dimensional quantum antiferromagnet with 5=1 from a 
spin-wave theory is Ah = 0.3914J (Ref. [27I ) which well 
agrees with the numerical density-matrix renormalization 
group (0.4105 J from Ref. l28h . The spin- wave theory for 
the one-dimensional quantum antiferromagnet with in- 
tegral spin can reproduce the static magnetic suscepti- 
bility which agrees well with the quantum Monte Carlo 
simulations. 2 - These previous results, together with our 
study, indicate that the formulation of the spin-wave the- 
ory is a valid formalism to apply to the low-dimensional 
quantum magnets. 

In summary, we have studied the QHAFM on a square 
lattice by using a perturbation spin-wave theory. The 
Lagrangian constraint for the absence of the local mag- 
netization plays an important role and leads to an ef- 
fective mass for the spin wave excitations at finite tem- 
peratures. In linearized spin-wave theory, the calcula- 
tion of the uniform static magnetic susceptibility can re- 
produce almost all features obtained from the quantum 
Monte Carlo and the high-temperature series expansions. 
It has a linear-temperature dependence at low tempera- 
tures and a broad peak feature at an intermediate char- 
acteristic temperature T c cx JS(S + 1). At high tem- 
peratures, the uniform magnetic susceptibility follows a 
Curie- Weiss law. The static magnetic susceptibility from 
our perturbation theory with the first corrections from 
the spin-wave interactions agrees well with the numeri- 
cal simulations. 
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